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A two-temperature Ising-Schelling model is introduced and studied for describing human 
segregation. The self-organized Ising model with Glauber kinetics simulated by Miiller 
et al. exhibits a phase transition between segregated and mixed phases mimicking the 
change of tolerance (local temperature) of individuals. The effect of external noise is 
considered here as a second temperature added to the decision of individuals who consider 
change of accommodation. A numerical evidence is presented for a discontinuous phase 
transition of the magnetization. 
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1. Introduction 

The problem of human segregation is an important problem of society and politics 
even in the 21st century Social sciences have been investigating the reasons 
and and nature of segregation for a long time. Sociologist have introduced several 
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models, one of them is the Schelling model . From physicist point of view that 
model is a 3-state voter-type non-equilibrium model (groups A,B and empty), with 
spin-exchange dynamics at zero temperature (T = 0) on a 2-dimensional square 
lattice. Although the model describes a segregation by a quench without external 
reasons, unwanted frozen states may also occur. 

Recently it was shown by computer simulations ^ that there exists a simpler 
model, namely the Glauber-Ising model which captures the essence of human 
segregation. Besides that, the usage of the simple T > spin-flip dynamics makes it 
computationally easier and one can avoid the frozen states of the Schelling model. 
In this model the temperature plays the role of a global tolerance; by varying it, 
the model may or may not evolve into the ordered (segregated) state. 

However a constant, global tolerance is rather artificial in a society, it can vary 
from individual to individual and can change in time as well. By introducing a local, 
time dependent temperature (tolerance), with a feedback mechanism from the local 
neighborhood the model becomes more realistic. The results do not change too 



1 



2, 2008 2:47 WSPC/INSTRUCTION FILE odor 



2 Geza Odor 

much a self-organization of the average temperature occurs. The parameters of 
this model are the global rate of forgetting (of tolerance) and the response of local 
tolerances on the neighborhood. The sum of these local changes determine the local 
temperature. Hereafter this self-organized segregation model will be called SO-Seg 
model. 

Stepping further towards more realistic models, one can pose the question what 
happens to this model if the decision of individuals are affected by an independent 
external noise as well. The external noise can be an artifact of a random environ- 
ment, housing, moving situation, presence of shopping centers ... etc. The external 
noise introduced here as a second temperature, i.e. individuals are connected to a 
second heath bath. 

Two-tcmperature two-state voter-type models have been investigated intensively 
recently and have become the prototypes of non-equilibrium models (for a review 
see ^ . An important finding of these studies was the discovery of relevant factors 
affecting the phase transitions of models exhibiting Z2 (up-down) symmetry ISl. In 
particular models with general, isotropic spin-flip dynamics maintaining the Z2 
symmetry can be classified as two temperature models, where one temperature 
controls the bulk, the other the interface fluctuations. The Ising model is a special 
case of these models, where both temperatures are nonzero and the time-reversal 
symmetry drives the system into an equilibrium state. The transition of these Z2 
symmetric, two-temperature models has been found to be continuous, Ising type 
unless the bulk temperature is zero. In the latter case it is first order, voter model 
class type In this work I investigate the effect of a second temperature applied 
as an external, independent heat-bath to the spins of the SO-Seg model. 

2. The model 

The model is defined on 2-dimensional square lattice, with periodic boundary condi- 
tions and Ising spins {si = (1,-1)) distributed initially randomly (zero initial mag- 
netization = no segregation). The kinetics follows a Glauber spin- flip (sequential) 
update with Heath Bath acceptance rate (see Ref.E]), depending on the local temper- 
atures. In the SO-Seg model each individual has four interacting nearest neighbors 
and a randomized initial local temperature, with an average value < Ti(0) >= 1.5. 
This local temperature is lowered by 5Tp at each update for modeling the loss of 
tolerance. This alone would just make a quench to Ti — > with domain coarsening. 

To model people's awareness of the dangers of segregation they can increase their 
own temperature (tolerance) by STa-, if all four neighbors of an individual belong to 
the same group as Si. If all four neighbors belong to the different group the local 
temperature is decreased by the same amount. 

Our external noise is described by a heat bath of a second temperature T2 
applied to the decisions. The actual spin-flip will be the logical OR of internal and 
external flip decisions. The magnetization M{t) =< Si{i) >, the average number 
of like neighbors minus unlike neighbors < N{t) >, and the average self-organizing 
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internal temperature (tolerance) < Ti(t) > is followed up to tmax — 2x 10^ Monte 
Carlo sweeps (MCS) of the lattice. 

3. Simulation results 

The simulations were performed on L = 400, 2000 and 4000 sized square lattices, 
up to tmax = 2 X 10^ MCS, with cooling rates: 5Tp = 0.01, 0.02 and tolerance steps: 
5Ta = 0.003, 0.002. The phase transition of the SO-Seg model is at 6Ta = 0.0029151, 
so 5Ta ~ 0.003 corresponds to super-critical, 6Ta — 0.002 to sub-critical situations. 

As Figure ([T]) shows the inclusion of a small second temperature T2 <^ 1 does 
not change the composition of neighbors and Ti in the steady state (as t 00). The 
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Fig. 1. Clusters survive small external noise: T2 = 0.5. 



same can be seen by plotting N{t) (see Fig. ([2])) and the average internal temperature 
of individuals (Fig. ([3])). For stronger external noise the domains are destroyed, but 
the average tolerance goes to zero too. This means that the unsegregated state 
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Fig. 2. Average composition as the function of time and external temperature T2 = 0, 1,2. ..10 
(top to bottom) 
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Fig. 3. The average tolerance of individuals as the function of time and external temperature 
T2 = 0, 1, 2... 10 (top to bottom) 



can be maintained with the help of strong external noise, without worrying about 
people's local tolerance. The same analysis for STa = 0.002 resulted in similar trends 
in N(t = 30000) (see Figg]) and in the average tolerance Ti{t = 30000). Note that 
one observe even a weak increase in the asymptotic values for T2 < 1. Running the 
simulations on larger sizes there were no change in this, excluding the possibility of 
finite size errors. 

By increasing the second temperature T2 the transitions of A^(oo) (see Fig. ([4])) 
and Ti(oo) are very smeared. The magnetization density m — Af(oo)/i^ on 
the other hand shows a sharp fall, indicating a first order phase transition at 
T2* = 0.135(2). A fitting attempt using the form m = A{Ti - T2Y did not re- 
sult in agreement with the 2d Ising class continuous phase transition behavior, 
characterized by /3 — 1/8 Furthermore a hysteresis cycle can also be found by 
starting the simulations with different initial conditions (ordered vs. disordered), 
which is a clear hallmark of a first order phase transition. This is in contrast with 
the results for two-temperature Z2 symmetric models exhibiting Ising transition in 
2d (the bulk noise in nonzero of course). One may understand the discontinuous 
transition here by realizing that this model is effectively a coupled system: Ising -f 
temperate field. The temperatures can increase if at least four neighbors are in the 
same state and decrease without condition. In the language of reaction diffusion 
systems it is a quadruple model, where at least for neighbors are needed for a cre- 
ation but the removal is spontaneous. It is well known that in 2d such a quadruple 
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model exhibits a first order transition . 

Simulations with different initial temperatures (ri(0) = 2.5), cooling rates {Tp = 
0.2) and 5Ta have not resulted in changes in the transition point T2*. 
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Fig. 4. Average composition as the function of external temperature {T2). Bullets (higher data 
on the left) correspond to 5Ta = 0.003 forgetting rate, crosses to 5Ta = 0.002 forgetting rate. 
Inset: The magnetization density for STa = 0.003 shows a first order transition. 

4. Conclusion 

A two-temperature, self-organized Ising-Schelling model has been introduced and 
investigated by numerical simulations. A low second temperature, which represents 
the external noise to the decision of individuals for moving does not change the 
segregation behavior of the model. A temperature bigger than T2 = 0.135(2) ran- 
domizes the segregated, ordered domains and results in low average tolerance of 
individuals. While the self-organized tolerance and the average composition of the 
steady state show continuous variation on T2 the magnetization exhibits a first or- 
der phase transition. The threshold did not show considerable dependence on the 
model parameters. 
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